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Abstract Strong gravitational lenses are unique cosmological probes. These pro- 
duce multiple images of a single source. Whether a single galaxy, a group or a clus- 
ter, extracting cosmologically relevant information requires an accurate modeling of 
the lens mass distribution. A variety of models are available to this purpose, never- 
theless discrimination between them as primarely relied on the quality of fit without 
accounting for the size of the prior model parameter space. This is a problem of 
model selection that we address in the Bayesian statistics framework by evaluating 
Bayes' factors. Using simple test cases, we show that the assumption of more com- 
plicate lens models may not be justified given the level of accuracy of the available 
data. 

Images produced by strong gravitational lenses result of different light-paths. 
If the source behind the lens has a variable luminosity, this will manifest with a 
time delay between the images. This time delay At depends on the gravitational 
potential of the lens, and the underlying cosmological model. Therefore, we can 
derive constraints on cosmological parameters (in particular Hq), provided a lens 
model is assumed. Hence, lens modeling as well as accurate measurements capable 
of discriminating between models are critical to the study of time delays. 

We aim to tackle this problem from the point of view of Bayesian model selection 
analysis (see e.g. iH). A large number of lens models have been proposed in a 
vast literature. Given the fact that observables are limited to the position of the 
images, their time delay and flux ratio, we restrict our analysis to simple examples 
characterized by a few parameters. In particular we consider two models for lenses 
with two images, so called "double" lenses (for a review on lensing, see [l]). 

1. Power-law model; assume a density profile p °^ r^", with n a free parameter. For 
« = 2, it describes an isothermal lens. In order to assess the dependence on the 



Irene Balmes 

Laboratoire Univers et Theories (LUTh), UMR 8102 CNRS, Observatoire de Paris, Universite 

Paris Diderot, 5 place Jules Janssen, 92190 Meudon, France 

irene.balmes@obspm.fr 



1 



Irene Balmes 




Prior on 7: <0.1<0.2<0.5 <0.'<0.2<0.5 <0.1<0.2<0.5 <0.1<0.2<0.5 

Prior on n; arqc smol orqc smol 

At only At and fijx ratios 

Fig. 1 Bayes' factor between model 1 and 2, with different priors. Above the dot-dashed line, the 
evidence in favor of model 1 is strong. Each color represents a different lens 



prior parameter interval we assume two different priors: < « < 3 (large) and 
1 < « < 3 (small). 

2. Power-law model with external shear: assume the previous model with the ad- 
dition of shear accounting for environmental effect on the lens. This adds two 
parameters: the strength of the shear 7, and its direction. Expected values for the 
shear vary up to 7 ~ 0.1, therefore we assume three different priors on 7: 7 < 0.1, 
< 0.2 and < 0.5 respectively. This allow us to test the shear strength up to nearly 
unrealistic values. 

We performed a UkeUhood data analysis for a sample of lenses and inferred the 
Bayes' factor for model 1 and 2 under different priors. 

Results are summarized in Fig.[T] Large Bayes' factors favor the simpler model, 
model 1. Above a certain threshold (dot-dashed line), the evidence in favor of 
model 1 is considered strong. In the following, we highlight a few relevant aspects. 

• Effect of the prior on «: The lens data set is mainly composed of galaxies, which 
we expect to be nearly isothermal. Nevertheless, our analysis shows that a large 
fraction of our sample is accurately described by model 1 if < « < 1 . 

• Effect of the prior on 7: In more than half of the cases, allowing higher (unreal- 
istic) shear strength does not change the Bayes' factor This illustrates the effect 
of the Occam's razor term in the Bayes' factor: a wider range for a parameter is 
bound to give a better fit, but this is balanced against a penalty factor. 

• Effect of the flux ratios: Time delays depend on the gravitational potential of the 
lens, whereas flux ratios depend on its second derivative. Furthermore, they are 
subject to a number of local phenomena that do not affect time delays. As a result, 
flux ratios require more complex models than time delays. This is consistent with 
our findings in Fig.[I] indeed, adding flux ratios as a constraint leads to having 
less lenses accurately described by model 2. 
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